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A return probability of random walks is one of the interesting subjects. As it is well 
known, the return probability strongly depends on the structure of the space where the 
random waker moves. On the other hand, the return probability of quantum walks, 
which are quantum models corresponding to random walks, has also been investigated 
to some extend lately. In this paper, we present a limit of the return probability for a 
discrete-time 3-state quantum walk on a hexagonal lattice. 
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1 Introduction 

Quantum walks are quantum analogies of random walks. As a popular topic in science, 
random walks have been intensively investigated. Quantum walks started to get attention 
in early 2000’s because their application to quantum computers was focused on around that 
time. We can consider quantum walks as quantum algorithms and applied them to some 
algorithms [1]. When we solve problems by using classically stochastic algorithms, some 
of them require exponential number steps for the size of their input to get the solutions. 
Since quantum algorithms give quadratically faster speed-up than classical ones, we expect to 
solve such difficult problems within polynomial number steps by using quantum algorithms. 
Quantum walkers on a lattice have coin-states and the behavior of the walkers is described 
by wave function. The wave function of a discrete-time quantum walk evolutes according to a 
unitary operator which is obtained from a product of a coin-flip operator and a position-shift 
operator. 

In this paper we treat a quantum walk on a hexagonal lattice. Hexagonal lattices show 
up in various fields of physics. For example, graphene has a hexagonal structure. To use 
quantum walks for applications, it would be important to clarify the behavior of quantum 
walks on such a realistic lattice. In particular, we concentrate on computing a limit of a 
return probability of the quantum walk after long time. It is well known that the return 
probability of random walks strongly depends on the structure of a lattice where the random 
walker moves. On the other hand, that of the quantum walker depends on its coin-flip 
operator rather than the structure. Our walker moves on a 2-dimensional rectangle lattice. 
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but the motion is controlled by position-shift operators and it is equivalent to a motion on a 
hexagonal lattice. There are a lot of limit theorems for quantum walks on a line, while a few 
2-dimensional walks have been exactly analyzed That’s because it is generally much 

harder to compute 2-dimensional walks than 1-dimensional walks. Kollar et al. [5] treated 
a 3-state quantum walk on a 2-dimensional plane. The walker defined by them moves on a 
hexagon-like lattice and they numerically analyzed the return probability. The decay order 
of the return probability of some special 2-dimensional walks were analytically estimated and 
discussed with the quantum Polya number [Sll]. Though it is a special class of 1-dimensional 
walks, a relationship between localization of the quantum walks and recurrence properties of 
the corresponding random walks was clarified [5]. 

This paper is organized as follows. In Sec. [21 we define a 3-state quantum walk on a 
2-dimensional lattice. The walker starts from the origin and its motion is controlled by two 
position-shift operators so that it moves on a hexagonal lattice. The return probability is 
analyzed and its long-time limit is calculated in Sec. El In the final section we summarize and 
discuss our result. Then we will briefly show a limit theorem in a rescaled space by time. 


2 Description of a quantum walk on a hexagonal lattice 


First, we define a quantum walk on a two-dimensional lattice H = Hi U 1HI2 with Hi = 
: a;, 2 / G z| ,H 2 = ^ where Z means the set of integers. 

Then we describe the motion of the walker on a hexagonal lattice by position-shift operators. 
Figure [H which will show up later, helps us get a picture of both the 2-dimensional lattice 
H and the hexagonal lattice which we treat in this paper. The position of the walker is 
expressed on two Hilbert spaces Kp-^ = {\x, v) : (x, v) G Hi} and = {|Xi it) : (x, v) G H 2 }, 
and a Hilbert space Hp = Hp^ 0 Hp^ just corresponds to the two-dimensional position space 
where the walker moves. In addition, the walker at each vertex on T-Lp supposes to be in 
superposition with three coin-states. We introduce a Hilbert space T-Lc which is spanned by 
the basis {|0), |1), |2)} so that the coin space of the walker is expressed. To compute a limit 
law later, we take the following orthonormal vectors: 


|0) = 

■f 

0 

, ll) = 

o' 

I 

, |2) = 

O' 

0 


0 


0 
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The whole state j^'*) of the quantum walker at time t G {0,1,2,...} is described on the tensor 
Hilbert space fip ® 'He- The position of the walker is shifted by two position-shift operators 
Si, S 2 after the superposition is operated by a coin-flip operator C as follows: 


|^i+i) = (5i + 52)C|^'t) 


( 2 ) 
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where 
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with c = COS0, s = sin0 {0 G [0, 27r)). Since the behavior of the walker is obvious, we will not 

i q = 2^ 

3’ ■’ 3 


treat 0 = 0,7r. When we set c = —s = the coin-flip operator C becomes a Grover coin 


C'= E ® 
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Assuming ('I'ol^'o) = 1, the probability that the walker is observed at position (x,u) G H, is 
defined by 


'[{Xt,Yt) = (x,u)] = {^t\ ( \x,v){x,v\^^\j){j\ ) 

j=o 


( 7 ) 


where (At, Yt) denotes the position of the walker at time t. Finally we set an initial condition 

|d/o) = |O,O)0(a|O) + /3|l) + 7|2)), (8) 
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for G C such that |ap + |/3p + I 7 P = 1, where C means the set of complex numbers. 

This means the walker starts from the origin, and if it takes to be the initial state Eq. ([ 8 |), 
the walk on a hexagonal lattice is obviously realized by the position-shift operators Si, S 2 - 
The position-shift operator S 2 (resp. Si) does not work to any whole system such 
that I'l') G Hpi (resp. G Up^), and since Cld^) G 'Hp^ (resp. G 'Hp,^), we see {Si -f 
5 ' 2 )C|'I') = 5'iC|'I') (resp. = 52(7 |d>)). Combining this fact and the initial state Eq. (| 8 ]), we 
should note that the whole system becomes |'I' 2 t) = {S 2 CS 1 C)* |'I>o) G Hp-^ and |d> 2 t-i-i) = 
SiC {S 2 CS 1 C)* I'f'o) G T-Lp^- The walker, therefore, is not observed on Hi at time 2t 1. 
Figure [T] explains the lattice H = Hi UH 2 and the vertices where the walker arrives according 
to the position-shift operators 5i, ^ 2 . 



from the origin O (0, 0), it moves over the points with a black circle. 

Moreover, as we can see from Fig. [2J the probability distribution sharply depends on the 
initial condition. 



(a) a = ^ = 7 = ^ (b) a = 7 = 0, /3 = 1 

Fig. 2. Probability distribution at time t = 100 (c = — 


3 Limit law of the return probability 

In this section we show a long-time limit law of the return probability. The return probability 
is studied as one of the interesting subjects in the field of classical random walks because it 
is involved with the recurrence. Similarly, the return probability of quantum walks is also 
interesting with the recurrence, compared with the classical ones. Regarding quantum walks 
on the 2 -dimensional square lattice, the relations between the return probability and the 
recurrence were discussed for some special walks mm- As we already defined the model in 
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the previous section, the walker starts from the origin. So, the return probability means the 
probability that the walker can be observed at the origin (i.e. P[(Xt,Yt) = (0,0)]). Since 
|^/>2t+i(0,0)) = ^[0,0,0] (i.e. P[(X2t+i,y2t+i) = (0,0)] = 0) because of | 'i'2t+i(o,o)) e "Hp^, 
we should focus on the probability at the origin at time 2t, where T means the transposed 
operator. For the return probability, one can get the following limit. 


Theorem 1 If the walker starts from the origin, we have a long-time limit of the return 
probability 


limP[(X2t,r2t) = (0,0)] = 

t^OO 
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Figure [3] shows how the return probability f‘[{Xt,Yt) = (0,0)] numerically converges to the 
limit Eq. 



time t 

(a) a = /3 = 7 = 

Fig. 3. The probability at 



(b) a = 7 = 0, /3 = 1 

origin at even time (c = — ^, s = ) 


Proof. We define a transform |'kt(a, b)) G (a, b G [—tt, tt)) of the walk at time t as 




^—i(ax+by) 
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V't (f, 4^) ) + E IV't (^, ) , (10) 
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If the walker starts from the origin, the transform is reduced to a simpler form, according to 
even or odd time: 
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That’s because we have |V' 2 t(Xi^)) = ^[0,0,0] at position (x,u) G 112 and |^/j 2 t+i(x,i^)) = 
"^[0,0, 0] at position (%, v) G Hi, for any t G {0,1,2,...}. The amplitude at the position where 
the walker can be observed, is extracted from a transform 



As we can see from this inverse transform, the amplitude is obtained just like the inverse 
Fourier transform. We, hence, call the transform |'I't(a, 6)) the Fourier transform in this 
paper. Equation ([2]) leads us to the time-evolution of the Fourier transform 


|'I'2t+i(a, b))=R{a,b)C\^ 2 t{a,b)), (15) 

|^ 2 t+ 2 (a, b)) =R{-a, -b)C\i/ 2 t+iia, b )), (16) 


where the matrix i?(a, 6) is a 3 x 3 diagonal unitary matrix 
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(17) 


(18) 


From Eqs. m and (HH), the Eourier transform at time 2t turns out to be 
|^ 2 t(a, b)) = {i?(-a, -b)CRia, 6)c}‘ |^'o(a, b)). 


(19) 


We express the eigenvalues Xj{a, b) {j = 1, 2, 3) of the unitary matrix R{—a, —b)CR{a, h)C as 
follows: 

A,(a,6) = e*‘'^(“-'') (j = 1,2,3), 


with 


i/i(a, 6) = 0, 

V 2 {a,b) = arccos {c^ — ^(1 — c)^ sin^ 6-|-cos a cos 6} , 
^^{a.b) = 27r — arccos {c^ — 1(1 — c)^ sin^ 6-I-cosacos6} . 

Here, for x, y G Z, we put 

i{ax+by) 


( 20 ) 


( 21 ) 
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2s2(l — cos a cos 6) -|- (1 — c)^ sin^ b 


( 22 ) 
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By using the residue theorem with regard to the parameter a in Eq. (I22L we get an integral 
representation 


gix,y) = / 

^0 


cos(&| 2 /|) |2s^ + (1 — c)^ sin^ 6 — (1 — c) sin6A/(3'+lip'^^7i'^^^0^'cos^| 
7r(l — c) sin b{ 2 s^ cos &)l“l ^(3 + c)^ — (1 — c)^ cos^ b 


-db. 

(23) 

Again, for x,y G Z, one can obtain a long-time asymptotic behavior of the amplitude at 


position ^ 





3x 
2 ’ 


2 


)) 


2 

(ni(a, b)\Ma, &)) ki(a, &)) ^ oo) 


'-§ {VEi(a, 7 )G(x, 2 /, 1,-1) + W2{a,l3)G{x -h 1, y - 1,1,-1) -h 112 ( 7 ,/3)G(a;, y-h 2,-1,1)}' 

= -^(1 - c) {lEi(a, 7 )G(x - 1, 2 / + 1,0,2) + lE 2 (a, /3)G(x, y, 0,2) + 1 ^ 2 ( 7 , /3)G(x, y, 0, -2)} , 

§ {lEi(a, -f)G{x, y, 1,1) -h W 2 (a, /3)G(x -11, y - 1,1,1) -h 1 ^ 2 ( 7 , /3)G{x, y, -1, -1)} 

(24) 

where |w_, (a, 6)) (j = 1,2,3) are normalized eigenvectors corresponding to the eigenvalues 
Xj{a, b) and we have put 


G{x, y, xi, yi) =g{x, y) - g{x - Xi,y - yi), 

11"i(z 1,-22) = - S2l -I- SZ2, 

\/2 

1E2(zi,Z2) =SZi -^(1 - C)Z2. 


(25) 

(26) 

(27) 


The asymptotic symbol hi(t) ~ /i2(t) {t —>■ 00 ) means limt_>oo hi{t) / h2(t) = 1. The Riemann- 
Lebesgue lemma have been used in Eq. (1241) . Straightforwardly computing the long-time 
asymptotic behavior of the amplitude at the origin, we have 


IV'2t(o,o)) 


'[ 1 _ V2sA(e) o f (3+c)A(6) il ■ 

1^2 / 7r(l-c) \ 7r(l-c) 2j' 

- {s« - \/2(l - c)/3 -F 57 } 

f (3+c)A(e) _ V2sA{e) o I 7i _ 

\ 7I-(l-c) 2/“ 7I-(l-c) \2 -K )l 

This also gives the long-time limit of the retnrn probability. □. 


(t —>■ 00 ). 


(28) 


4 Discussion and summary 

The return probability of symmetric simple random walks on a lattice strongly depends on 
the structure of the lattice, while that of quantum walks depends on their coin-flip operator. 
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It have not been cleared yet how the return probability of quantum walkers is determined 
by the structure of a lattice. We treated a quantum walk on a hexagonal lattice and its 
coin-flip operator was given so that it included a Grover coin. Since the walker started from 
the origin, the return probability was simply defined by the probability that the walker can 
be observed at the origin. As a result, we derived a long-time limit theorem about the return 
probability. Also, Theorem [T] gives a condition that either localization or delocalization at 
the origin occurs. When we get limsupj_^oo P[(Ar 4 , Yj) = (x,u)] = 0 (resp. > 0) for position 
(x, y) £ m, let us say that delocalization (resp. localization) at position (x, v) occurs. From 
this definition, we find a condition of delocalization at the origin. Delocalization at the origin 
is realized if and only if 


, , o V2(l + c) 

|a| = —-—, fi = - a, 7 = a. 


(29) 


The condition of Figs. [5]- (a) and[31-(a) satisfies Eq. (1^ under the condition c = — i,s = 

On the other hand, a limit theorem on a rescaled space by time t can be also demonstrated. 
Theorem 2 There exists a continuous function f : i —> K. such that, for a:, y £ R, we 

have 

2Yt 


lim 
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<y = 


r ,u f 

^ — oo J — ( 


dv X{a,j3,^)5o{u,v) + f{u,v), (30) 


where 5o{x,y) is the Dirac 6-function at the origin and 
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— 1 > 5ft(a7). 


(31) 


We mean that M is the set of real numbers and 5i(z) {z G C) is the real part of the complex 
number z. 


We can also consider the value A{a,f3,^) as a return probability on the rescaled space. 
Although this theorem is obtained by using the Fourier analysis like the other past researches 
(e.g. liiaii]), it is not completed because the continuous function f{x, y) is not computed. 
Note that we have A(q:,/ 3,7) = 0 for the initial state such that Eq. (l2^ . It would be one of 
the interesting future problems to compute the continuous part f{x, y). 
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